Dynamic correlations in an ordered c(2x2) lattice gas 
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We obtain the dynamic correlation function of two-dimensional lattice gas with nearest-neighbor 
repulsion in ordered c(2x2) phase (antiferromagnetic ordering) under the condition of low concen- 
tration of structural defects. It is shown that displacements of defects of the ordered state are 
responsible for the particle number fluctuations in the probe area. The corresponding set of kinetic 
equations is derived and solved in linear approximation on the defect concentration. Three types of 
strongly correlated complex jumps are considered and their contribution to fluctuations is analysed. 
These are jumps of excess particles, vacancies and flip-flop jumps. The kinetic approach is more 
general than the one based on diffusion-like equations used in our previous papers. Thus, it becomes 
possible to adequately describe correlations of fluctuations at small times, where our previous theory 
fails to give correct results. Our new analytical results for fluctuations of particle number in the 
probe area agree well with those obtained by Monte Carlo simulations. 
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INTRODUCTION 



Lattice gas models have been extensively studied (see 
[l|). They are used for the description of various physical 
phenomena such as surface mass transport, ionic con- 
ductivity, etc. and are equivalent to the familiar Ising 
model with conserved spin dynamics (e.g. Kawasaki dy- 
namics Q)- The general theory of kinetic phenomena 
in lattice systems is far from complete, despite its long 
history, started from seminal papers 0, @]- Additional 
complications and new effects arise in particular models 
due to nontrivial lattices and complex interactions. All 
this explains the permanent activity in the field ranging 
from keystone problems such as the susceptibility of two- 
dimensional Ising model to more specific questions 
like those considered in our previous studies j5|, i6j, 0] and 
other recent papers [1, Q . 

The subject of the present paper is kinetic phenomena 
in two-dimensional lattice gas with nearest neighbor re- 
pulsion in the ordered c(2x2) phase (see Fig. [1]), which 
corresponds to an Ising antiferromagnet. Static proper- 
ties of this system are well studied by means of appro- 
priate approximate or exact methods. Exact expressions 
for the spatial multipoint correlation functions in grand 
canonical ensemble at half- filling [10] arc the most pro- 
nounced theoretical results obtained in this field. In con- 
trast, temporal correlation functions arc less studied even 
in the case of thermodynamic equilibrium (excluding the 
simplest cases like on-site correlation function ^ ) . While 
most efforts were undertaken to describe critical dynam- 
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ics, still kinetic phenomena outside the critical region are 
poorly understood. The standard approach for classical 
gas based on BBGKY hierarchy is not directly applica- 
ble to the case of lattice gas, because the equation for 
one-particle distribution function involves not only two- 
particle but also higher distribution functions. 

To describe kinetic phenomena, we use the method 
of essential configurations which can be outlined as fol- 
lows. There are powerful cluster expansion methods like 
virial or low-temperature expansions for studying static 
properties of the ordered equilibrium phase. The main 
idea of these methods is to reduce the total configura- 
tion space to some important (essential) configurations 
which involve only finite number of structural defects. 
To study the dynamics of lattice gas, some additional 
configurations, which are generated by already chosen 
ones in course of their evolution, should also be consid- 
ered. These transient configurations, having short living 
times, do not contribute essentially to static properties 
but affect the kinetics because of their participation in 
the displacements of structural defects. Thus, by taking 
into account mutual transformations of both defect types 
we may derive kinetic equations in the reduced space of 
essential configurations. Our main interest is for the two- 
point two-time correlation function. It will be obtained 
from the solutions of the mentioned kinetic equations. 
This function is a building block for macroscopic time 
correlation functions which play a central role in kinetic 
phenomena 11]. 

From the point of view of practical application the 
calculated correlation function can be utilized in theory 
of adsorbates [l^, [13] • Modern experimental techniques 
give the opportunity for detailed investigations of the mo- 
tion of such small objects as individual atoms on crystal 
surfaces [l3] (it is worthwhile to mention that not only 
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kinetics but even femtosecond dynamics can be studied 
[lH). Thus, calculated correlation functions of particle 
density fluctuations in scales of one-particle jump may be 
compared with those experimentally obtained. In such 
a way some microscopic parameters of the adsorbate- 
adsorbate and adsorbate-substrate interactions may be 
obtained from this comparison. In fact, our paper is an 
attempt to generalize the so-called fluctuative method 
to the case of kinetic scales. The correlation func- 
tion is also connected directly with the dynamic structure 
factor measurable with diffraction techniques [l^,[l^ . Fi- 
nally, the theoretical analysis may be used to understand 
and to interpret data obtained by means of Monte Carlo 
simulations. 

More specifically, the objectives of the paper are as 
follows: (i) to develop an analytic method for studying 
kinetic phenomena in ordered lattice systems; (if) to sup- 
plement the approach used in the recent paper [5] ; (iii) to 
explain previous results ^] of MC simulations for short- 
time correlations of fluctuations in small probe areas. 

The paper is organized as follows. In Section |TT] we 
specify the system under consideration and define the 
main quantities. Section IIIII serves as general outlook 
to the system: typical snapshots, structural defects, 
and some insight into kinetic phenomena. Then, in 
Section IIVI the description of essential configurations, 
derivation and solutions of the kinetic equations for struc- 
tural defects will be presented. In Section |V] the correla- 
tion function is calculated. And finally, the comparison 
of analytical data with the results of MC simulations is 
given in Section fVll 



II. SYSTEM SPECIFICATION 

Consider a lattice gas of N particles hopping on a 2D 
square lattice of size Lq x Lq with periodic boundary 
conditions (i.e., on a torus). Points of the lattice (sites) 
are denoted by a single letter, e.g., x — {xi,X2), where 
xi^2 — 0, . . . , Lq ^ 1. Double occupancy of a site is for- 
bidden. The Hamiltonian is 

where Uxy = C/i if a; and y are nearest neighbors and 
zero otherwise. The case of repulsion is considered {Ui > 
0). We denote pUi = 0, where /3 = T^Ms the inverse 
temperature. The system is described by two parameters: 
the average concentration of particles c (in surface science 
it is usually denoted by 9) and the interaction strength 
0. 

Also, it is assumed that: 

(i) the system is at subcritical temperature, i.e., ip > 
(j>c ~ 1.76 so that q = e~'^ is a small parameter; 

(ii) c is nearly 0.5 so that c(2x2) phase is pronounced; 

(iii) the system is at thermodynamic equilibrium 
(which is to imply the absence of domain walls). 



Under these conditions the system is represented by 
two sublattices (Fig. [T]): one is almost empty and the 
other is almost filled. They are distinguished by indices 
a, /9, . . . and a, r, . . . for the empty and filled sublattices, 
respectively. The concentrations of excess particles in the 
empty sublattice and vacancies in the filled one are small 
enough to treat the system as rarefied gas of structural 
defects. 

The motion of particles is assumed to be realized by in- 
stantaneous jumps. Thus, simultaneous displacements of 
different particles are forbidden (single-particle-jump ap- 
proximation) . The hopping itself is an overcoming of an 
activation barrier. It should be noted that in real physi- 
cal systems as in surface adsor bates, all microscopic pa- 
rameters including adsorption sites, hopping pathways, 
activation barriers, adatoms interaction, etc. can be cal- 
culated by ab initio computations (see e.g. [20j). 

To be specific, we assume the rate of individual parti- 
cle jumps from a filled site x to any unoccupied nearest 
neighbor site y to be given by W^y = I'oe'"^ , where k is the 
number of filled sites, which are nearest-neighbors to the 
site cc, and Vq is the jump rate in the absence of particle- 
particle interaction. The detailed balance conditions are 
satisfied for this choice of jump rate: Wxy/Wyx = e'^'""'^", 
where Sx is the particle energy in site x expressed in units 
of T. Hence in the absence of external perturbations, the 
system relaxation towards equilibrium is ensured. 

This "single site energy" probability, which is widely 
used in modern surface science, has been suggested in 
pioneer paper [2l| , in which the results of computer sim- 
ulations of particle migration for two dimentional lattice 
gas were reported. The same model was used later to 
develop analytical approaches in [H, [2^ . It is believed 
now to be a more realistic representation of the diffusive 
kinetics in experimental systems, although it is slower 
than the METROPOLIS algorithm. 

Without the loss of generality we choose the unit of 
time in such a way that vq = 1. Thus the motion of a 
single particle in the empty lattice are described by the 
following master equation: 

Px^ -^Px + ^Px+e, (1) 
e 

where the sum runs over nearest neighbor sites (e means a 
unit vector). Here and throughout the paper, the length 
is given in units of a lattice constant. 

The main object of the paper is the equilibrium two- 
point two-time correlation function given by 

{5nx{t)6ny{Q)) = {nx{t)ny{Q)) - {nx{t)) {ny{Q)) . (2) 

In practical applications other quantity is used: the cor- 
relation function of particle number fluctuations in a 
small probe area, {SN{t)SN{0)) . It can be measured ex- 
perimentally and useful for surface diagnostics. These 
two correlation functions are connected by the relation 

{dN{t)dN{0)) = {Snx{t)6ny{0)) , (3) 
x,y 



where the summation is over the probe area. It is more 
convenient to use "per site" quantity, which in case of a 
square probe area with L x L size is defined as 
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{SN{t)SN{0)) 

l2 



(4) 



In the hmit L — > cx) (in such a way that Lq ^ L) we 
obtain the thermodynamic correlation function S(t) = 
Soo{t), which for t = can be calculated also by the 
formula 5(0) — T (9c/9/i)y where /i is the chemical 
potential. 

The dynamic structure factor can be also expressed via 
two-point correlation function: 

SUt, k) = j^Yl {SnAt)Sny{0)) e-'''^^-y\ (5) 

so that 5^(^,0)= SL{t) (this normalization differs from 
that used in p^.[l9l|). 



III. GENERAL OUTLOOK 

A typical snapshot (obtained in the course of MC sim- 
ulations) is shown in Fig. [TJ The "chessboard" phase is 
clearly seen, even for cf) close to its critical value. The 
"defects only" view is presented in Fig. [21 where only ex- 
cess particles in the empty sublattice and vacancies in the 
filled one are shown. This view supports the idea that 
we can treat the system as rarefied gas of structural de- 
fects. To elucidate the nature of these defects we redraw 
the snapshots, leaving only the topological charge. This 
quantity is defined as the difference between the number 
of excess particles and vacancies belonging to the same 
cluster of connected defects in Fig. [21 To produce Fig. [3] 
we define a cluster as a group of defects which are nearest 
neighbors, and place the topological charges at the mean 
coordinates of the clusters. The result is shown in Fig. [3l 
from where it is clearly seen that the seeming variety of 
the structural defects in Fig. [2] actually represents the 
fluctuating deformations of some basic defects of charges 
1, —1, and 0. Also sparsely populated pairs of defects 
(the large circle in Fig. [3|) are present. 

Before classifying these basic structural defects, it 
should be noted that in a real system domain walls are 
often presented as nonequilibrium long-living metastable 
configurations. Nuclei of such inclusions are always 
present as large continuous areas of defects, observed in 
Fig. [21 The true domain walls make the whole system 
essentially inhomogeneous, lacking good ergodic proper- 
ties. Therefore, we assume that there are no domain walls 
in the neighborhood of the probe area. 

Various structural defects may be classified by their 
concentrations at exact half-filling. To get rid of possible 
ambiguity in defects definitions, we use the following rule: 
two structurally similar defects are called distinguishable 
at level q" if their concentrations at exact half-filling dif- 
fers by 0{q^), k < n. To make the ground for the study 
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FIG. 1: A snapshot for (j> = 2.0 and c = 0.5. A black square 
is an excess particle (occupied site), while a white square is a 
vacancy (empty site). 
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FIG. 2: "Defects only" view of Fig. [T] The full circles rep- 
resent excess particles, while the empty circles are vacancies. 
The sites which are filled according to the regular "chess- 
board" pattern are not shown. 
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FIG. 3: Snapshot of the topological charges after the subtrac- 
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of kinetic phenomena we will consider the lowest-order 
defects shown in Fig. 2] up to level q'^. It is worthwhile 
to note that in Fig. [2] (where q = 0-8qc) only five clusters 
have smaller concentrations: q'^, q'^, q^, q"^ , q^^. 
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FIG. 4; All structural defects up to level: (a) excess parti- 
cle, (b) flip-flop pair, (c) excess particle monomer, (d) flip-flop 
tetrad, (e,f) transient configurations during monomer jump, 
(g) isolated flip-flop pair. Only essential sites are shown, po- 
sitions with respect to sublattices can be easily guessed. Con- 
figurations obtained by color inversion must be also included. 

First of all, excess particles and vacancies can be easily 
distinguished as alone circles in Fig. [51 The concentration 
of both defect types at exact half-filling is given by [l| 



1 - 



16g2 



(1-9)^ 



g2 + 4g3 (6) 



(for simple cubic lattice n = z; « Q)- Only these 
two structural defects constitute the level because all 
other defects either have concentration o{q^) or are in- 
distinguishable from these two at the level q^. 




FIG. 5: Structural defects at q'^ level, from left to right: excess 
particle monomer, flip-flop pair, vacancy monomer. Essential 
sites are encircled. 

At the level q'^ of classification we distinguish excess 
particle surrounded by only occupied sites and excess 
particle with vacancies in the nearest neighborhood (see 
Fig. [5]) . In the first case the structural defect is called ex- 
cess particle monomer, Fig.Ul^c). In the second case only 
one vacancy is allowed to have concentration distinguish- 
able from zero at level g"^. This is so called "flip-fiop" 
pair, Fig. [lUb), that is a pair of adjacent excess particle 
and vacancy. Their concentration per site of a sublattice 
is 



2 (nove) « 2(7^ 



(7) 



Note that this value is considerably greater than in the 



no 



= {no) - {noVe) 



q'-5q\ 



(8) 

The difference between excess particle, Fig.SJJa), and its 
monomer. Fig. [DJc) can be illustrated by the equation 



from which also the notion of "essential sites" used in 
the figures can be understood. Only monomers and "fiip- 
fiop" pairs constitute the level g'^. 

The structural defects of the level g"* are monomers, 
isolated fiip-flop pairs, flip-fiop tetrads, and transient 
configurations arising in the course of monomer jumps 
(Fig-IHc-g), see also Fig.[TT|). Despite the low concentra- 
tion of transient configurations q^), they are innegli- 
gible for description of monomer motion. Starting from 
this level one must consider also monomer pairs, that is 
two-particle configurations. At g^-levcl this reduces to 
taking into account onsite generation-recombination of 
vacancy and excess particle monomers and onsite exclu- 
sion principle for identical monomers. 
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FIG. 6: Generation-recombination precursors. The structural 
defect (a) is also depicted in the box to show its position with 
respect to sublattices. 

The generation-recombination processes must be con- 
sidered for understanding the long-time asymptotics of 
the correlation function (IjSBl)- The generation passes 
in three steps. The first two are the creation of a double 
fiip-fiop pair resulting in one of the transient configura- 
tions shown in Fig. [6] (their concentration is ~ g^). At the 
last step adjacent vacancy and excess particle monomers 
are created. 
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FIG. 7: (a) Isolated monomer, (b) side and (c) corner dimers 
of excess particle. 

At the level q^ isolated monomers become distin- 
guishable from side and corner dimers shown in Fig. [T] 
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FIG. 8: Relative contribution of long-distance components of 
the correlation function to S{0) at exact half-filling accord- 
ing to Eq. |9] depending on interaction strength q/qc = e'^""'^. 
Parameters of MC simulations (see Sec. IVIfl are marked by 
crosses on the horizontal axis <j) ~ {3.0,2.7,2.4,2.0,1.86} 
(from left to right). 

The concentrations of these structural defects have lead- 
ing terms of the low-temperature expansion — ISq*, 
q'^—Sq^, and q'^ — 9q^, respectively. At this level the short- 
range attraction between monomers arises, that can be 
observed as small excess in concentration of dimers com- 
pared to squared concentration of monomers, which is 
^ q'^ — IQq^. The number of dimer configurations in- 
creases when c ^ 0.5 (as squared monomer concentra- 
tion) and in that case they may contribute significantly 
to mass transport [5i, |6|, lD] . 

Obviously, when approaching the critical point or mov- 
ing away from exact half-filling the larger structural de- 
fects should be accounted to obtain reliable results. To 
understand the importance of large defects for lattice gas 
statistics we calculate their contribution to S'(O), using 
the analytic expression for (SrixSno) — ((5nj;(0)(5no(0)) at 
exact half filling (the expression is not explicit but can be 
derived by recurrent procedure [Kl]). This can be done 
by evaluating 

^"(0)= {^^i^.^.)^^m) (9) 

\xi-\-X2 I >77l 

and comparing it with S{0) = S'*'(0). The result is shown 
in Fig. [5] for m — 3 and 9. Roughly speaking, the upper 
curve, in which short-distance correlations (|a;i-f a;2| < 2) 
are excluded, gives an estimate of an error arising when 
only monomers and fiip-flop pairs are taken into account. 
This means that the contribution of long-distance cor- 
relations in the case of large values of (j> is negligible. 
Thus, our approach seems quite reliable when > 2.5 
(T < 0.7 Tc). 

For further analysis it is useful to distinguish two time 
scales of the kinetic phenomena. At small time scale we 
have a diversity of structural defects (typical snapshot is 
shown in Fig.[T]). The major part of them have very short 



lifetime. Their fast evolution is due to strong repulsive 
interaction with the nearest neighbors. The contribu- 
tion of short-term fluctuations to the dynamic correlation 
function decays exponentially with time. At long times 
only topologically stable structural defects (monomers of 
excess particles and vacancies, as well as domain walls) 
survive and contribute to the correlation function. The 
motion of long-living defects results in such macroscopic 
phenomena as diffusion, segregation, domain growth. A 
typical decay of the corresponding correlation function 
is given by a power-law form. In the next section we 
will derive kinetic equations governing the evolution of 
structural defects. 



IV. KINETIC EQUATION FOR STRUCTURAL 
DEFECTS 

The configuration space of the system is extremely 
large, 2 ». But observations made in the previous section 
together with general principles of statistical physics sug- 
gest that only a tiny part of the entire configuration space 
is actually occupied. In the present paper we derive the 
kinetic equation in the reduced space of essential config- 
urations. The general sketch of the method is as follows. 
At first we set the desired concentration level, say q'^, and 
choose configurations with concentration no less than q'^ 
(note that each configuration has Lq translationally in- 
variant copies, therefore we use "per site" quantities). 
Thus, we obtain the configuration space of structural de- 
fects. Then we determine transient configurations which 
may be considered as "deformations" of the structural 
defects in the course of their displacements. It should 
be emphasized that the search of the transient states is 
not a specific problem of lattice systems exclusively. It is 
a general problem of kinetics. And finally, we write the 
standard master equation for the Markov chain in the 
reduced space of essential configurations (defects -I- their 
deformations). We have to solve this set of equations 
up to order g*^. If we are interested in the correlation 
function we may exclude transient configurations (with 
concentration o{q^)) from this set. The equations ob- 
tained in such a way govern the evolution of structural 
defects only. This evolution is nonmarkovian and this is 
the main effect of transient states. 

In our case the first nonvanishing level is q^, which 
includes monomers only. We will also consider the next 
level, q^ , that is, fiip-flop pairs in order to look for the 
corrections due to higher orders defects. 

The motion of monomers was explained in 6]. The 
main idea of that paper may be outlined as follows. Let 
us consider an excess atom shown in Fig. [S] In such 
configuration the most frequent processes are jumps like 
that from Di to an empty site with rate e"^. Let it will 
be site C. The new configuration is short-lived because 
the most probable process has rate e"^"^. Now there are 
two possibilities. The first is that the displaced particle 
returns to site Di so that the final state coincides with 
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the initial one (i.e. Di — > C ^ Di). The second is that 
the particle in site A jumps to site Di and the final state 
is one with the excess particle displaced from A to C 
(i.e. Di ^ C, A ^ Di). Both possibilities have equal 
probabilities. Therefore if one neglects the living time of 
the transient state, the rate of this complex (two-step) 
jump of the defect will be given by ^e"^ [6|]. Similarly 
the rate for A — > B complex jump is 2 x ■^e'^, where the 
factor of two arises from two possible transient states: Di 
and D2. The same reasoning for vacancy displacements 
(Fig. [T0|) results in rate value equal to 1 for A ^ B jump 
and i for A — > C jump. 
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FIG. 11: Essential configurations for excess particle monomer: 
(a) excess particle monomer, (b,c) transient states. Centering 
is indicated by arrow. 
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FIG. 9: Displacement of excess particle monomer. 
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FIG. 10: Displacement of vacancy monomer. 



The essential configurations at level include: empty 
c(2x2) lattice, excess particle monomer (see Fig. [TIJa)), 
vacancy monomer, and transient states for both types of 
monomers (see Fig. [TlTb.c) for excess particle). Other 
transient states may be ignored. Because the empty lat- 
tice and the two kinds of monomers have different topo- 
logical charges, monomers move independently at this 
level of accuracy. Let us now derive the kinetic equa- 
tion for excess particle monomers. The approach may 
be generalized to any dimension. Here we consider d- 
dimensional case that makes formulae more obvious. The 
level should be replaced by q"^. 

Let us denote the probability of the configuration in 
Fig. Illf a) by Pa , where a is a site indicated by the ar- 
row in the figure. The transient states in Fig. [TlTb.c) 
are denoted by p^*^ . They are formed by vacancy at 
site (7 and two excess particles at sites a + e and a + e' 
fFigs. fnTbl andfTlTc) correspond to e' 7^ e and e' = — e, 
respectively). Now the kinetic equation for excess parti- 



Ee 



Pa+e'), e'^e. 



cle monomer is given by 

fp„ = -2d(2rf-l)g-Va + 
\ pf = -2gi-2'ip-' + q-\p,+e + 

(10) 

In the first equation of this set the first rhs-term describes 
(a)^(b,c) processes. Here 2d is the number of outer par- 
ticles in monomer and 2d — 1 is the number of sites to 
which they can jump with the rate q~^. The second rhs- 
term corresponds to the reverse processes which all have 
the same rate q^~^'^. In the second equation of the set the 
first rhs-term describes (b)— >(a) (if e-|-e' ^ 0) or (c)— >(a) 
(if e + e' = 0) processes. Here the factor of 2 accounts 
for two possible particles which can jump to the site a. 
The last term corresponds to the reverse processes. 

To obtain the correlation function we have to solve this 
set of equations with respect to Pa up to the order 0{q''). 
To do this we introduce a new variable 



Po 



p;~ 



(11) 



and also a deformed Laplace operator (in continuous limit 
it reduces to standard laplacian) 



4d 



x-\-e-\-e' 



(12) 



Now (|T0|) reduces to two compact equations: 



-2d(2d-l)g-Va + gi-2''K, 



Ad{2d - l)g"Va + ■idq^^Aap. 

(13) 

They can be solved by means of discrete d-dimensional 
Fourier transform: fx f{k) = ^}Zx f^^^^^ ■ inverse 
is given by 



fx 



(27r)« 



/(fc)e 



dfc. 



(14) 
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In the fc-domain (|13p transforms into two ordinary hnear 
differential equations of the first order: 



p = -2d{2d - l)q-^p 

f>' = -2gi-2V + Ad{2d - l)q-^p - 4:q-^6{k)p, 



where 



d{k) = d^ 




(15) 
(16) 



Their solution is expressed via two exponents e^^-^*, 
where A1.2 are solutions of the equation 

A2 + 2A [d{2d - 1)9-1 + qi-^rf] ^ 45-2<i^(fc) = 0. (17) 

At small q these solutions are given by Ai w —2q^^'^'^, 
describing nonpropagating short-time fluctuations, and 
A2 ~ —2q~^5{k), corresponding to random displacements 
of a monomer as a whole. 

Similarly, we may perform a Laplace transform of (I13|) 
with respect to variable t. In such a way Green's func- 
tion of Eq. [12] corresponding to pa (that means pa {t) = 
J2f3Ppi^)G'a-l3i*)) reduces to 



sq 



2d 



Ad 



sq 



1 



sq 



2d 



Ad 



+ d{2d-l)q 



2d-2 



(18) 



where s is the Laplace variable, tilde marks transformed 
function, upper index "e" means "the excess particle 
monomer" , and g denotes Green's function of the equa- 
tion Lfx = Aa;(p, which in fc-domain reduces to 



g{t,k) 



exp^i 




(19) 



Simple estimate of shows that G^_^(i) — 

ga^-a {2dtq-^) + 0(g2d-2). 

The kinetic equation for vacancy monomers is similar 
to that for excess particles but with the only difference 
that time t is rescaled by the factor q. This means that 
q~^ and (^~'^'^ in should be changed by 1 and q^~^'^ , 
respectively. Therefore Gl_^{t) = .g^_r (2dt) + ©(g^-i-S). 
The diffusion coefhcients of monomers, derived from the 
mean square displacement — Afc(7(i, A:)|fe=o — 2dDt, are 
equal to 2dq^^ and 2d for excess particles and vacancies, 
respectively. 

Summarizing, in the two-dimensional case at level q'^ 
the evolution of monomers is governed by the following 
Green's functions: 

G^-/5W=.9a-;3(4ig-i), G^_, (i) - .g,-. (4t) , (20) 

where g is given by its Fourier transform (|19p with d = 2. 
Note that gx is nonzero only if |a;i -I- . . . -I- x^l is even. 



This simply means that monomers move only on their 
own sublattices. 

For |a— /3| ^ 1 and |(t— r| 3> 1 we come to the diffusion 
equation used in our previous studies [1, H, 0] • 



O 
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FIG. 12: Essential configurations for fiip-fiop pairs. Centering 
is indicated by arrow. 

Essential configurations for flip-flop pairs include only 
empty c(2x2) lattice, see Fig.[T2ja), and flip-flop pair, see 
Fig. HHb) . We denote the probability of these configu- 
rations as Pq and p^, respectively. The kinetic equations 
are simple: 



Po = -2dpo + 2dq'-'"'pe, 

Pe = -gl-^'^Pe +P0- 



(21) 



Here — 2c?po term describes (a)— >(b) processes, where 2c? 
is the number of nearest neighbor sites to which the par- 
ticle may jump. 

The corresponding Green's functions are given by 



"-"00 — 



2d + i^ 



Gff 
Oe — 



I _ ^-(2d+u)t 

2d + v ' 



rff _ „rff rff _ 1 2d + z.e-(^'^+'-)^ 2d-l 



Gff 



2d 2d+v 
1 2d-|-i^e-(2''+'')* 1 
2d^ 



2d 



2d + v 



2d 



e + e', (22) 



where v = q^"^"^. It follows from Eq. [2T]that in the case 
of g <g; 1 the relaxation time of the fiip-fiop fiuctuations 
is given by g^'^"^. 



V. CORRELATION FUNCTION 

Now we can obtain the correlation function. There is 
no correlation in the motion of the two types of monomers 
and fiip-fiop pairs. They evolve independently and give 
additive contributions to the total correlation function. 
To find these contributions we can use the general for- 
mula for the correlation function of independent random 
walks, derived below. 

Let N particles randomly walk on a lattice. Denote by 
the position of i-th particle. The number of particles 
in site x is given by 



N 



(23) 



where is the indicator of an event A. In what fol- 

lows we take into account the identity of particles. In 
this case their average number at a given site, that is the 
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average on possible trajectories with fixed initial distri- 
bution, reduces to 

M)) = N {I{at) = x}) = NP{at) = x}^ Np^it), 

(24) 

where is the position of any selected particle and the 
formula {T{A}) — P{A} is used. The correlation func- 
tion can be calculated as follows {t > s): 

i 

+ 5](ZU.(t)=x,0(s)=y}) 

i^j (25) 
= 7V(XU(<)-x,e(s)=y}> 

+ N{N - 1) {x{m = ^}> inm = y}) 

= Npy^is, t) + N{N - l)p,{t)py{s). 
Hence the correlation function will be given by 

{6n^{t)Sny(s)) = {n.j;{t)ny{s)) - {n^{t)) {ny{s)) 

= Npy{s)[GyS~s)~p,{t)]. (26) 

For homogeneous lattice at equilibrium reduces to 

{5n.,{t)5ny{s)) = n[G,_,,(t - s) - G,_y(oo)], (27) 

where n = {nx{oo)) is the equilibrium concentration of 
particles. 

Straightforward application of (I27p to monomers gives 
their contribution to the correlation hmction to be 



(28) 



{5n„{t)5nr{Q)Y = n^Gl_,{t) = nV_^(4i), 
where the concentrations of monomers were calculated in 



(29) 




r ± c - o 



(upper sign is for excess particle monomers). Other com- 
ponents like (dnaSua-) are zero for monomers in the ap- 
proximation explained above. 

For flip-flop pairs the factor n in P7|) is unity (it is 
the concentration of conflgurations in Fig. fTWa)'). We 
should take also into account that each particle in the 
almost empty sublattice may appear there due to a jump 
from any of its nearest neighbor sites. In this way we ob- 
tain the following nonzero components of the correlation 
function: 



(30) 



VI. FLUCTUATIONS IN PROBE AREA: 
COMPARISON WITH MC SIMULATIONS 

Now according to ^ we have to sum up the derived 
two-point correlation functions over the square probe 



area of size L x L to obtain the quantity SL{t), which 
we can compare with the results of MC simulations. 

From the contribution of monomers will be given 
by 



STit) 



L2 



^5a-,3(4ig-i) + ^^ff.-.(4t). (31) 



a,/3 



The sum in JSTj) can be evaluated by using the identity 

2 



■^'^ \ — -rr 'rr]d 



dfc, 



(32) 



which takes place if x and y vary within the same domain. 
In particular, for parallelepiped (xi = 0, Li — 1) 



-ikx 



sm 



kiLi 



n-sini 

2—1 2 



(33) 



In case of even L the sums in (f3T|) and ((32|) are connected 
by the identity X;„,/3 = | T,x.,y ... 

For flip-flop pairs we can obtain explicit expression by 
using the formula (time arguments are omitted) 



'(00)''«(00) 



) 



L-1 



'(mO)f"-(00)y 



m—l 
L-1 



+ E ~ '"^^ ('5"(mm)'5"(00)) 



(34) 



+ 4EE(^--)(^-0('5r 



m— 1 

L-1 m-1 



L2 

m=2 1=1 

Thus, from ^U^) it follows 



'(m^)''f^(oo)/ 



(35) 



This expression does not depend on c at level q'^. The 
total correlation function is given by 



Sr — S'f 



(36) 



To analyse the contribution of flip-flop jumps to fluctu- 
ations it should be noted that Sf is a decreasing function 
of size L. This point may be easily understood by the 
observation that displacements of only border particles 
(whose total number is proportional to L) determine fluc- 
tuations of SN. Relative contribution of flip-flop jumps 
to S'l is shown in Fig. 1131 from which it follows that this 
contribution is negligibly small when (/) > 2.5 (see also 
Fig. [H]) . This result agrees well with the calculations 
shown in Fig. [S] Having zero topological charge, flip-flop 
pairs contribute insigniflcantly to the correlation func- 
tion at times of the order or longer than q^ . On the other 
hand, Fig. [Ol shows the adequacy and good accuracy of 
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FIG. 13: Relative contribution ol fluctuations due to flip-flop 
jumps for (j} = 2.7 (top) and 3.0 (bottom), L = 6. Circles 
and diamonds are MC data, while full lines are theoretical 
calculations. 

our approach. The main shortcoming of that be- 
comes apparent for small (j) in Fig. ll3[ originates from the 
crude underestimation of concentration of flip-flop pairs 
in Q. 

For large values of interaction parameter (0 > 2.5) 
the motion of monomers gives dominant contribution. 
It was shown in [5] that in the case of large times the 
generation-recombination processes of monomers should 
also be taken into account. These processes are responsi- 
ble for establishing local equilibrium in the defect system 
and determine the dissipation of smooth spatio-temporal 
inhomogeneities of the particle density and their fluctu- 
ations. Bringing together the approach of the present 
paper, which provides adequate description of the sys- 
tem at kinetic scales (the characteristic length is of the 
order of lattice constant and the characteristic time is 
of the order of defect living-time in a given site), with 
that developed in 5] for hydrodynamic scales we get a 
description of fluctuations at any time, both short and 
long. Thus, it becomes possible to compare theoretical 
data with the results of MC simulations in the whole 
range of computer simulations, Fig. [W) A good agree- 
ment between the two can be seen. 



VII. DISCUSSIONS 

A very simple physical picture explaining the nature 
of dynamic correlations in the ordered c(2x2) phase may 
be outlined from this analysis. In the case of suffi- 
ciently strong particle-particle interaction (0 > 2.5 or 
T < O.TTc) and half-flUing (|c - 1/2| < q) the dom- 
inant contribution to the dynamic correlation function 
is due to random walks of two types of structural de- 
fects: excess particle monomers and vacancy monomers, 
which have concentrations q^. At the same time, an- 
other two processes must be accounted for accurate de- 
scription of the correlation function: 1) at short times 



'L 
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0.001- 
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FIG. 14: Correlation function of fluctuations in probe area 
per site vs. MC steps for (j> = 2.70, c — 0.5, L is indicated. 
1 MCS corresponds to time [4(1 + e'^)(l + e^'^)]-\ Circles 
are MC data. The solid line at short times is the correlation 
function of fluctuations due to random walk of monomers only 
given by Eq. 1311 The dash line (merging with solid at 20 
MCs) is the same quantity with flip-flop jumps accounted as 
well fEg. I36p . The solid line at long times is the correlation 
function of fluctuations from the diffusion and recombination 
of monomers from Ref. jS]. 



- additional fluctuations caused by flip-flop pairs with 
concentration 2q'^; 2) at long times - faster decay of cor- 
relations due to generation-recombination of monomers, 
O(g^). If |c— 1/2| ^ q then one of two types of monomers 
is in majority with concentration |2c — 1|, and its ran- 
dom walking determines the correlation function. In this 
case dimers with concentration |2c — must also be 
accounted, which decrease the collective diffusion coeffi- 
cient. 

The decay of correlations is described by several char- 
acteristic times. First of all, it is the living time of flip- 
flop pair (Fig. [TWb)). which is of the order of q'^. Two 
other characteristic times are connected with the dura- 
tion of monomer jumps (g^ for excess particles and q^ 
for vacancies). They are equal to living times of tran- 
sient configurations (Fig. Illf b.c)) formed in the course 
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of two-step defect displacements. The relative contri- 
bution of these transient configurations is of the order 
of q'*. At this time scale all other short- living config- 
urations snapshoted in Fig. [5] also decay. At moderate 
times only monomers (and their groups) are essential. 
The corresponding characteristic times are q for excess 
particle monomers and 1 for vacancy monomers, which 
are time intervals between their successive jumps. Ran- 
dom walk of monomers results in power-law decay of cor- 
relations, as t ^ GO. At very long times generation- 
recombination processes become essential. They speed 
up the decay of correlations though the asymptotic be- 
havior of the correlation function is still t~^. This is in 
full accordance with some rigorous lower estimates (Q 
and ref. therein) and with long-time asymptotics of re- 
laxation of concentration fluctuations m A + B C re- 
versible diffusion-limited reaction considered in (see 
also [H, [26] for details). 

Summarizing, the proposed method of essential con- 
figurations makes it possible to describe the evolution of 
structural defects and to obtain the dynamic correlation 
function for a lattice gas with nearest neighbor repulsion 
in the ordered c(2x2) phase. Our calculations explain 
MC simulations reported in Q at short times and show 
the range of times where the diffusional approach devel- 
oped in [1, 0, 0] is not applicable. 



Further development of the present work includes the 
extension of the method to level, that is necessary 
for accurate description of generation-recombination pro- 
cesses. The main complication is that we must proceed 
from one-particle description to a many-particle one. 

Concluding, our consideration shows that for known 
particle-particle interaction and jump rate mechanisms, 
the correlation function of particle number fluctuations in 
small probe area can be calculated for the specific model 
of the ordered lattice gas. Hence the comparison of the 
analytical and experimental data concerning short time 
correlations becomes possible. In principle, the informa- 
tion about individual particle jumps or complex jumps 
accompanied with the defect displacements may be ex- 
tracted from such comparison. This paper may be con- 
sidered as the attempt to get better understanding of 
short-time correlations in ordered systems. 
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